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A calculation of the ^He transverse (e, e') inclusive response function, Rt, which includes A de- 
grees of freedom is performed using the Lorentz integral transform method. The resulting coupled 
equations are treated in impulse approximation, where the NNN and NNA channels are solved 
separately. As NN and NNN potentials we use the Argonne V18 and UrbanalX models respectively. 
Electromagnetic currents include the A-isobar currents, one-body N-currents with relativistic cor- 



■ rections and two-body currents consistent with the Argonne V18 potential. Rt is calculated for 

the breakup threshold region at momentum transfers near 900 MeV/c. Our results are similar to 
those of Deltuva et al. in that large A-isobar current contributions are found. However we find that 
these are largely canceled by the relativistic contribution from the one-body N-currents. Finally a 
comparison is made between theoretical results and experimental data. 

PACS numbers: 25.30.Fj, 21.45.-v, 14.20. Gk, 21.30.-x 
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I. INTRODUCTION 



It is well known |l| that subnuclear degrees of freedom play an important role in nuclear 
dynamics. In conventional low-energy nuclear physics the relevant subnuclear degrees of 
freedom are considered to be mesons and nucleon isobars. Electron scattering affords an 
excellent tool for studying these degrees of freedom which are manifested in the transverse 
response through meson exchange (MEC) and isobar (IC) currents. The consideration of 
such subnuclear currents has a long history in the physics of few-nucleon systems. For 
two-body systems a review can be found in j^. One important issue in the MEC is their 
consistency with the nucleon-nucleon (NN) potential being used. Such consistent MEC's 
have not only been taken into account in deuteron electrodisintegration, but also in the 
electrodisintegration of three-body systems js-?]. On the other hand IC have not received 
the attention in three-body systems which they have in the two-nucleon sector. Nevertheless 
there exists a rather complete calculation by js], Isl Iq] where and A degrees of freedom 
have been treated on an equivalent level via a coupled channel calculation with NNN and 
NNA channels. Also A-effects have been studied (sl in ^He electrodisintegration below the 
three-body breakup threshold using the transition-correlation-operator method jlO |. 

The present work incorporates the dynamics of the A-resonance into the many-body 
wavefunction by means of the impulse approximation (lA) [11]. This method, as in the 
transition-correlation-operator and coupled channel methods, avoids the static approxima- 
tion by fully including the kinetic energy in the A-propagator. For electromagnetic deuteron 
breakup it has been shown that the A-effects resulting from an lA calculation are rather 
similar to those resulting from a coupled channel calculation if the energy is sufficiently 
below the resonance position of the A {12I. 

A calculation of Rt in our case requires an integration over continuum states of the 



coupled NNN -|- NNA system. As has been demonstrated previously [13|, [IJ] the Lorentz 
integral transform (LIT) method is well suited for calculating inclusive quantities such as 
response functions. Examples of its use in calculating electron scattering response functions 
of three- and four-body nuclei employing realistic nuclear forces (two- and three-body forces 
include (i) nonrelativistic calculations of Rl for three- nucleon systems Q, for ^He Q 
and inclusion of relativistic effects in and '^He [17| and (ii) calculations of Rt with 

n p n 

relativistic and consistent MEC contributions for three-body nuclei [6|, |7|, |18[ . There is an 
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LIT calculation [19] of Rt for ^He but with semi- realistic A^A^ forces. The method has not 
previously been applied to the coupled NNN+NNA system so that this paper is the first in 
that regard. 

The paper is organized as follows. In section II we describe the general formalism includ- 
ing the incorporation of A degrees of freedom in the LIT formalism. Section III specifies 
the input to the dynamical equations developed in section II. This includes subsection A 
detailing the potentials used in the NNN and NNA sectors, subsection B describing the elec- 
tromagnetic current operators used and subsection C outlining calculational details. Finally 
our results are discussed in section IV. There we compare our results for Rx to those from 
another calculation and to experimental data. 

II. FORMALISM 

In the one photon exchange approximation the cross section for the process of inclusive 
electron scattering on a nucleus is given by 



where Rl and Rt are the longitudinal and transverse response functions respectively, 00 is 
the electron energy loss, q is the magnitude of the electron momentum transfer, 6 is the 
electron scattering angle, and Q"^ = q"^ — u"^. 

In the present work we study the transverse response function. 



with A degrees of freedom within a non-relativistic approach. The low u region is considered. 
In ([2]) the subscripts i and / label, respectively, an initial state and final states, and the 
matrix elements are taken between internal states, center-of-mass motion being excluded. 

As mentioned in the introduction we employ the lA in order to take into account the 
A-resonance. This approximation is used for both the '^He ground state and the final state. 
Below we outline the various theoretical aspects required to include A degrees of freedom 
in a calculation of Rt via the LIT method. 

We consider the three-nucleon system with A^(939) and A(1232) degrees of freedom. 



d^l duj 




(1) 





3 



which leads to the following Hamiltonian 
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H = T + Sm + V = J2iTi + Smi) + J2^ij^ (3) 

i=l i<j 

where Tj is the kinetic energy of particle i with mass rrii, Snii = mi—m^ is its mass difference 
with nucleon A^(939), and Vij = Vji is the potential between particles i and j. By omitting 
the contribution from more than one A-isobar excitation, we construct the three-particle 
bound state from an NNN part and a NNA part i.e., 

l^o) = |^n + l^o^)- (4) 
The wave function is determined by the Schrodinger equation 

(Tn + V^^^ - ^)|^^) = -V^^^.^^ 1^^) , (5) 
iHA-E)\^^) = i6m + TA + V''''-E)\^^) = -V^^^'^^K) , (6) 

where ^^(^2,^1^2 = ^3^^. V,j{NiN2 ^ N[N^), = ^ yNA,NA^ 

^2^=1 '^i is denoted by T/v and Ta for the NNN and NNA channels, respectively. One 
should note that \7^^ is different from the usual sum of realistic A^A^ potentials, V^'^ , 
because the latter already contain implicit effects due to the A (e.g., in meson theoretical 
A^A^ potentials realized via part of the a meson exchange). 

Here we do not search for a direct solution of the coupled channel problem represented 
by Eqs. ( 15|6|) . Instead we use the lA where one computes and separately. More 
specifically one first determines the NNN part by solving 

{H^-E)\^^) = (7) 

with 

i7^ = T^ + \/^^ + \/^^^, (8) 

where V^^^^ = ^^^^-^^ Vijk is a three-nucleon force. In the lA one then uses the solution 
1^^) in order to calculate I^^q) through ([6]). 

Treatment of the continuum in the LIT technique requires the calculation of a localized 
Lorentz state |\^). This Lorentz state |\E') also has NNN and NNA parts written as 

1^) = |^^) + |^^). (9) 
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These fulfill the coupled equations 

(T^ + -Eo- = -y^^'^^ l^rA) + Oj,j,\^^) + O^vaI^o^) , (10) 

^Sm + T^ + V''^-Eo-a)\^^) = + OatvI^o^) + OaaIO , (11) 

where Eq is the three-body ground-state energy, the complex a = + iaj is the argument 
of the LIT in the transformed space, and the On^n2 denote the various diagonal (A^i = A^2)) 
and transition (A^^i 7^ A''2) electromagnetic current operators. One first solves for the NNN 
part using if at of Eq. (IHl) : 

{Hj,-E,-a)\^>'') = -K^^'^^(i/A-i5o-^)-MOAA.|^^) + OAA|^^)) 

+0jvivlO + 0jvAlO. (12) 

The above equation is derived by solving ( |TT|) formally for |\E'^), inserting the solution in 
(ITOl) . and dropping the term V^^'^^{H/^ — Eq — a)~^V'^^'^^ , since, as mentioned, A-effects 
to the nuclear interaction are already contained in the realistic H^q. With thus obtained 
one then calculates |^^^) in a second step through ffTTl) . Given the solutions |\^^) and |^^^) 
the LIT is obtained from the norm of the Lorentz state as 

(^1^) = (^^^1^^) + (13) 

These two terms correspond to different contributions. It can be shown that the piece 
(^^1^^) describes contributions due to final states with nucleons only. In this case the A 
degrees of freedom only contribute as virtual intermediate states. On the contrary the term 
^\[fA|v[fA^ describes contributions from final states containing a real A. The contribution to 
Rt from this term vanishes below the threshold for A production. Since the present study is 
for energies below that threshold this term will not contribute here to Rt- Such a real A has 
to decay into a nucleon and a pion eventually, thus the contribution corresponds 
to resonant pion production. 



III. INPUT TO DYNAMICAL EQUATIONS 



A. Potentials 



In the pure nucleonic sector we use the Argonne V18 (AV18) NN potential 20] and the 
Urbana IX (UIX) NNN potential [21|] while for the pure NNA sector we take ^^^^^=0 as 
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in the lA calculation of 12 1 . The NNN and NNA sectors are coupled via the V^^'-'^^ and 

n 

YNA,nn potentials. We use the same form for this potential as described in [12| except 
that here we use the short range cutoff given in the AV18 20| potential. In detail we take 
V{NN — )■ A^A) between particles 1 and 2 to have the form 



Vl2 — O^NN'n'T 1-T2 



Vo{m^) H -Vo{mp) 



O-NN'tt 



+ 



O'NN'iT 



S\2 



with 



3 

1 + — 

mr 



3 



mr 



:i-e 



(14) 

(15) 
(16) 



where cTj (r^) are regarded as transition operators for spin (isospin) of particle i, the coupling 



constants a^N'-w and a^N'o are taken from 



12|, c = 2.1 fm ^ is the same value as in the 



AV18 potential, and r = |ri — is proportional to the Jacobi vector (see Eq. (l22l 



B. Electromagnetic Currents 



In Eqs. (ITOl) and (ITT]) the driving terms are the transverse electromagnetic currents acting 
on the ground state. The term Oj^^ represents the purely nucleonic one- and two-body 
currents. For these the same one-body and two-body currents as employed in [3] are used. 
There the one-body currents included relativistic corrections to order and the two-body 
currents were consistent tt- and p-MEC currents constructed using the method of Arenhovel 



and Schwamb 



221 ]. The other terms Oat a, Oatv, and Oaa are transition and diagonal one- 



body A-isobar currents. The currents involving the A-resonance are given in Fig. 1. For 
one-body A-isobar currents we use the forms 

3 ./ AAT 



jAAf 



k=\ 



X qj 



2mA 



G 



AN AN 
Ml 'zk ) 



(17) 



and 



JAA 



where r', = 



k=l 

Rem, and p'^ 



2pl. + ^A 



2m/ 



Geo + 



i{ak^ X q) 
2mA 



G 



A 

Ml 



1 + r^k 



(18) 



Pfc — Pcm^k/A are the relative coordinate and momentum 



operator of the k-th particle, while Rem and Pcm are the center-of-mass coordinate and 



wwji 



VWVA/ 



J J 

NAA/W' WVW 



(a) (6) 
FIG. 1: Operators involving the A-resonance in the NNN (a) and NNA (b) channels 

initial total momentum variables of the system. With the assumption that Pcm is directed 
along q the term k in ( fT8|) has the value k = 1 + 2PcmAk/{Aq). However since here we 
are dealinsr with transverse currents this term does not contribute. The formfactors for the 
above currents are the same as in and take the form 

4.59 



'Ml (Q 



2x _ (4.35/2.0) 



(19) 



with A^ATi 
6 
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3m^ (1 + QV^Ini? 
840 MeV and AAAr2 = 1200 MeV . As in our previous NNN calculations 



18| we use the approximations 

G%AQ')^ut{QL)GW) 



l^p (Qav) 



av ) 



Ge^QIv) 

A/l( 



G^niQa 



GliQ 



2 ) 

avJ 



C. Calculational Details 

In order to solve Eqs. ([7]), (ITT]) and (fT2|) for the ground state and Lorentz vectors we 
expand the bound and Lorentz states on a complete antisymmetric basis. The reason for 
antisymmetrizing NNA states is that they couple to purely antisymmetric nucleonic states 
through symmetric operators. Thus the excitation of an antisymmetric NNN state to a NNA 
state occurs via an operator symmetric with respect to nucleons. Such an operator is a sum 
of operators which replace a nucleon with a A which therefore leads to an antisymmetric 
NNA state. For the NNN part we take the same correlated hyperspherical basis as in our 
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previous three-body calculations without considering A degrees of freedom (see, e.g., jl^). 
For the part with one A-excitation we use the following hyperspherical basis 

Wk) = ^^-^^=M= I [h^KiLDc ® {{siS2) Ss,) S)^j^^ , {{Uh) Th) TMt) , (20) 

where N is the order of the hyperradial function, K is the hyperspherical angular mo- 
mentum, L and / are the orbital angular momentum of the pair and of the spectator, 
respectively, coupled to the total orbital angular momentum C. Individual spin (isospin) 
quantum numbers of the three particles are denoted by Si {ti), i = 1,2,3 while the pair 
spin (isospin) is denoted by 5* (T), the total spin (isospin) by S (T) and Aij- stands 
for the projection of the total isospin. Quantum numbers and Aij denote the total 
angular momentum and its projection and (... (g) ...)jmj denotes {CS)J coupling. The 
index k denotes collectively {N, K, L,l, C, S,S, J', Aij,T,T, Air} B. We define the 
quantity i? to be if the pair of the three particle system contains one A and therefore 
(si, S2i S3) = (ti, ^3) = (1/2) 3/2, 1/2), and to be 1 if the pair contains no A and therefore 
(■51,52,53) = (^1,^2,^3) = (1/2,1/2,3/2). Note that if particles 1 and 2 are nucleons we 
always assume that L+S+T=odd so that the NN pair is already antisymmetric. The spa- 
tial basis functions in coordinate representation are products of hyperradial functions and 
hyperspherical harmonics 

{$,li2\hNK(Ll)CMc) = RN{p)YKCMc(^^^^^'P^^^^^'f^)^ (21) 

The coordinates p and x are defined in terms of the Jacobi vectors 



as p = -|- ^^^"^ ■, X = (^1 — ii)/ P- The coordinates Xi = cosOi and cpi are spheri- 
cal coordinates of the unit vectors in the directions of and ^2- We use the notation 
Ai = mi/rriN. Note that particle permutations entering the antisymmetrization operator 
interchange not only particle position vectors but also their mass numbers Ai. One has 
p"^ = + C,2 = ^1^1 + ^2^"! + ^s'^s " AR^, where R is the center-of-mass position. Thus p 
remains invariant under particle permutations. 

The operator (1 - (12)) / [2(1 + B^)]^^'^ makes the (12) -pair explicitly antisymmetric. 
Note that it gives unity if the pair contains no A but rather is an antisymmetric NN pair. 
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Finally the operator (1 + P)/-\/3, where P = (123) + (132), makes the three particle states 
with antisymmetric (12) pair totally antisymmetric. It turns out to be convenient to keep 
both B = and i? = 1 in (!20|) (thereby resulting in an overcomplete basis) and to fi- 
nally select out numerically the linearly independent states. This enables one to select out 
those states which give negligible contributions to the results. Application of the operator 
(1 + P)/V3 in (EO]) results in the more practical form 

\v>k) = L^lt2, (^^-^ ® ((^1^2) S's',) S)j^^ , {{t[Q T%) TMr) , (23) 

where j denotes collectively {B' S'T' , BSSTTN K{Ll)} and 



\FjCMc) — ~~7^^{^B',B^S',S^T',T\hNK{Ll)CMc) 



~^9b'S'T',BSSTT \ ^'nK{LI)CMc) ~^ 9b'S'T',BSSTT \ ^NK{Ll)CMc) ) (^^) 



with {Ll)C coupling for the total orbital angular momentum C and its projection Aic- As 
mentioned above components with B'=l in Eq. fl23l) represent configurations NNA in which 
particles 1 and 2 are nucleons. Those components with B'=0 represent NAN configurations 
in which particle 1 is a nucleon and particle 2 is a A. More details of the spin-isospin 
factors g'^ and g'^, and the spatial functions h'^ and h'^ are given in Appendix El Techniques 
employed in calculating the kinetic energy and the A^A^ — )■ A^A or A^A — )■ A^A^ potential are 
given in A|)pendices [B] and [C] respectively. 

As in [6i| all currents are expressed in terms of multipole expansions. Explicit expressions 
for the multipoles of the one-body current (containing relativistic corrections) are given 
in 18|. The multipoles for the tt- and p-MEC are found in [6] with modifications due to 
the implementation of consistent MECs for the AV18 potential listed in 3]- Finally the 
multipoles required here for the one-body currents relating to the A are listed in Appendix 
Id] With these multipoles one can then decompose the LIT of the response function according 
to its multipole content as 

^T(a) = -i^ J2 Y.^2J + l){RTyj\ (25) 

A=el,mag j] 

where 

(RtYj' = {^%\^%) = {^j'mI^j'm) + (^Jl/l^^iJ)- (26) 
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Here is the initial state total angular momentum, J and M. are the final state total 
angular momentum and its projection. The 1^;^) are the solutions of (11111121) where the 
following replacement is made on the rhs of these equations 

0|v^o) ^ W^m) = [T' ® l^o(Jo))]^A^. (27) 

In Eq. fl26l) A4 is arbitrary. In Eq. fl27p above \l/o is either or while O represents 
the various electromagnetic current operators on the rhs of fllO|lip . By projecting the rhs 
of (|27|) on the basis states (|23|1 one obtains 

/f^ I A A {(FfC'®{{s'A)S's',)S')j,\\T>^\\MJo)) 
{{Fj'c ® {{siS2)S s^)S)j,^,\q'j^) = Sj'^j6m',m -j^ ■ 

(28) 



We use the Lanczos method to calculate the response function, as described in 23|]. The 
response function is thus calculated by using 

- -Tm(0^^ I ^ \o^^ ) 

a J njQ -\- a — 

= _^l/^(gy^|0^)(0^| ^ |0n)(0n|Qy;vi), (29) 



m,n 



N 



where Q corresponds to the rhs of (fT2|) and |0„,) is the set of orthogonal Lanczos vectors. As 
starting vector \(f)Q) we choose the rhs of f|T2|) at one particular value of a. The expression 



{(pmliEo + a — HN)^^\(t)n) can be calculated by using Eq. (71) of [23|. A difference from 
previous LIT applications appears in the potential term on the rhs of the LIT equation, 
namely the A^A transition potential V^'^'^^ in (IT^ . The contribution of this term to 
{(l>n\QjM) is given by 

- Yl (0nl¥'^')A^z"fc'(¥'fc'|V^'^'^''^'^l¥'^)iV^^'(¥'^l0n^)(0m|(i^A ' Eo - a)-'\q^j^j^) , (30) 

I'k'lkm 

with 

\q^m) = \rL, ® \<{J^))]jM + [tL, ® \'^o{M]jM, (31) 

and the norm matrix A^ = {{LpkWi)). Also the second term in Eq. ( 12^ . can 
be calculated in a similar way with the Lanczos method. 

Because in this paper we are working at low energies near the ^He breakup threshold only 
the lowest multipole transitions contribute. We found sufficient accuracy by restricting the 
maximum value of to 5/2. The LIT is computed with o"/ = 5 MeV. The LIT inversion 
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24j | is made with our standard inversion method [ij, |25|. As discussed in p] we subtract 
from the LIT of the Ml transition the elastic contribution and invert the remaining inelastic 
piece separately from the other multipole contributions. 

IV. RESULTS 

In the present work we have used the LIT method to calculate A-IC effects on the 
transverse electron scattering response function RT{q,oj) for q^900 MeV/c and u up to 
20 MeV above the breakup threshold. This is the first application of the LIT method to 
include A degrees of freedom in the calculation of inclusive (e,e') response functions. The 

n 

importance of A-effects at these kinematics has previously been shown by [5|]. As A^A^ 
and NNN interactions we used the AV18 and the UIX potentials respectively. Following 
the lA calculation of [l^ we do not consider a diagonal NA interaction, i.e. V^va = 0. 
For the '^He ground state, our interaction model leads to a A-probability of 1.14 % which 
compares to 1.44% obtained by who used a CDBonn+A coupled channel potential model 



26|. In addition to the A-ICs jATA, Jan and Jaa the purely nucleonic currents include the 



nonrelativistic one-body current with relativistic corrections up to order [18] and an 

MEC consistent with the AV18 potential j?]. Concerning the relativistic corrections we leave 
out the u dependent relativistic piece in the present work as its contribution is negligible in 
the threshold region we consider. For the neutron magnetic and the proton form factors we 



take the dipole fit while the neutron electric form factor is taken from [27j. 

Fig. 2 displays our Rt results for several calculational options. The dominant transition 
multipolarity contributing at these near threshold energies is Ml. One sees that relativistic 
effects reduce the Ml transition strength considerably. If in addition MEC are also taken into 
account then the Ml contribution drops markedly leading to a rather different low-energy 
behavior of Rt- Inclusion of A-IC restores some of this lost strength and demonstrates, as 
anticipated, that the A effect is quite large. 



Now we turn to a comparison of our results with those of p]. For the comparison one 
should keep in mind that there are differences between the two calculations. Thus in [5]], (i) 
relativistic currents have not been considered, (ii) the Coulomb force is neglected in the final 
state interaction, (iii) their nucleonic potential model, the CDBonn 



is different from ours 



and does not reproduce the '^He binding energy and (iv) the full coupled channel calculation. 
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co.^^ [MeV] 



co.^^ [MeV] 



FIG. 2: Theoretical results for Rt of ^He as function of internal excitation energy u:\at at g=862 
and 927 MeV/c with various current operators: non-relativistic nucleon one-body (dotted), with 
additional relativistic corrections (dashed), plus additional MEG (dash-dotted), further addition 
of A-IC (solid). 



CDBonn+A [26|, is a more consistent treatment of A degrees of freedom than the lA, but 
leads to a slight underestimation of the '^He binding energy. There is another point which 
makes the comparison a bit more difficult. Namely, the Rt of [Sj is not calculated for a 
constant momentum transfer, in fact q is slightly decreasing with growing energy. Therefore, 
in Fig. 3, we prefer to display the results for each q in two panels. One sees that despite 
the various differences mentioned above the A-effects in both calculations are very similar. 
However, one also notes that relativistic corrections lead to an opposite effect, which is of 
the same size at g = 862 MeV/c, but somewhat weaker at g = 927 MeV/c. By comparison 
of results for Rt at about q = 500 MeV/c from [7] against those of sl one finds again an at 
least partial cancellation of relativistic and A contributions close to the breakup threshold. 
The stronger ..crease of ■„ the ealeulafo. of at h,gher energies, seen m Fig, 3. partly 
originates from the non-constant momentum transfers used in [5j as mentioned above (see 
also discussion of Fig. 4). 



Finally in Fig. 4 we compare the results of our calculation with experimental data 29|. 



Again, as in the calculation of [5j] the momentum transfer is only quasi-constant. Therefore, 
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FIG. 3: Rt of ^He. Lower panels {q = 862 MeV/c left, q = 927 MeV/c right): theoretical result 
from present work with non-relativistic nucleon one-body current and MEG (dotted) and additional 
A-IG dash-dotted, further addition of relativistic corrections for nucleon one-body current (solid). 
Upper panels {q at threshold as in lower panels, but slightly varying with uJint, see text) : theoretical 



results from Deltuva et al. 



5|] with non-relativistic nucleon one-body current and MEG (dotted) 



and additional A-contributions (dash-dotted). 

in the left panel we show Rt including all current contributions for the two extreme q values, 
i.e. q = 862 and 850 MeV/c, and in addition for q = 862 MeV/c the result where only A-IC 
are left out. The lower q corresponds to the data at about Umt = 20 MeV, while the higher q 
corresponds to the threshold energy. In the right panel of the figure we only show the results 
for q = 927 MeV/c that corresponds to the data close to threshold. The A-IC contribution 
is seen to be essential for obtaining a good agreement between theory and experiment below 
10 MeV. However, at higher energies the increase due to A-IC is not sufficient to describe 
the data even if one considers the slight shift of q with higher energies represented by the 
dotted curve. The present case is rather similar to deuteron electrodisintegration at higher 
momentum transfer, where at low excitation energy the leading Ml transition also has a 
minimum. For the deuteron case it is known (see e.g. 31|) that various theoretical ingre- 
dients, like for example the potential model dependence, can lead to rather large variations 
of the theoretical result. Thus our present study cannot give a final answer concerning the 
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comparison of theory and experiment. 




co^^ [MeV] 



co.^^ [MeV] 



FIG. 4: Rx of '^He for same kinematics and the same dash-dotted and solid cm'ves as in Fig. 2; 
in addition: result with all current contributions al q = 850 MeV/c (dash-double dotted curve). 



Experimental data with slightly varying g-values from 



We summarize our work as follows. We have illustrated how A degrees of freedom 
are integrated into the LIT formalism for a calculation of the inelastic inclusive trans- 
verse (e, e') response function Rt of '^He. The resulting coupled equations for the Lorentz 
states of the NNN and NNA channels contain, as opposed to the corresponding coupled 
Schrodinger equation, source terms with electromagnetic operators acting on the nuclear 
ground state. The A degrees of freedom are present in three different forms: (i) in the 
potentials V^^i ^2.^1^2^ ^^j^g A-propagator, and (iii) in the current operators jAr^Ar-^. The 

coupled channel equation is solved in impulse approximation, where the NNN and NNA 
channels are treated separately. First, the NNN part is solved using a realistic nuclear 
interaction with A^A^ and NNN potentials. The result thus obtained is then used for the 
solution of the NNA channel. The former gives a contribution to the electrodisintegra- 
tion of a purely nucleonic final state, whereas the latter leads to a contribution to the pion 
production channel. In the present work we have studied A-effects in Rt of ^He close to 
the breakup threshold at an momentum transfers of about 900 MeV/c. The response func- 
tion is affected by sizable MEC contributions, and, as in a previous full coupled channel 
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calculation [5| we find a considerable increase of Rt due to A degrees of freedom. Unlike 
the calculation of jsl we here take into account relativistic corrections to the nonrelativistic 
one-body current operator. At the kinematics considered here these relativistic corrections 
nearly cancel the A-IC contribution. This cancellation in fact leads to good agreement of our 
theoretical Rt with experimental data at very low energy transfer, while the experimental 
Rt is underestimated at somewhat higher energies. 
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Appendix A: Details of g'^ ,h'^ , g'^ ,h'^ 



The spin-isospin factors are 



9b'S'T'BSSTT 



V2 



T^f B' S'T' ,BSSTT 
B'S'T',BSSTT 





Si-H-S2,tl-(-S>t2 



9b'S'T'BSSTT 



1 \S'+T' 1 f 
— -Lj ■^JB'S'T',BSSTT 





if B' = 1 and B = 0, 
a B' = and B = 0, 

otherwise. 

a B' = 1 and B = 0, 
a B' = and B = 1, 
otherwise. 



fB'S'T',BSSTT = {-)^^'^S'Sf'f { ^. 1 ! i.„ )■ 6,'_,J,>^,J,L... 



s[s',S']^(t[t',T' 
The spatial functions in coordinate representation are 



(Al) 



(A2) 



(A3) 



NK{Ll)CMc 



) = (-) {^'li^l, ^2)^2(^1, ^2)\hNK{Ll)CMc)^ 



(A4) 



{^A2\hNK{Ll)CMj ~ i^'li~^1^^2)^2i~^1^^2)\hNK{Ll)CMc)^ (^5) 

where $,'i{^i,$2) ^2(^15^2) connected to and ^2 through the cycle operator (123) 
by {r[r'2r'^s[s'2s'^t[t'2t'^\ = (r2r3riS2S3Sii(:2i':3^i|, and using (!22|) we have 
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^'1(^1,^2) 



{A2 + As)iAi + A2) 

{A, + A2 + A3)A2 
{A2 + A3){Ai + A2) 



^1 



^1 



^2(^1 + ^2 + A3) 
{A2 + As){A^ + A2) 



(A2 + A3)(Ai + A2: 



^2- (A6) 



Appendix B: Kinetic Energy Calculational Details 



For a basis of antisymmetric states the kinetic energy can be written as 



A1 + A2 



-TT 
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(Bl) 



2mjv A1 + A2 + A3 

where tti = —id/d^i is the Jacobi momentum conjugate to $,1. Noting that for 
calculating the matrix elements of Tzf between the basis states ( p3|) one may drop 
(l-(12))/[2(l + 5"2)]V2^weget 



JV B"S"T" 1 2 I 3 

X {Fb'' S"T\B' S'ST'TN'K'(L'V)CMc\'^\Wb" S"T" ,BSSTTNK{Ll)CMc) ■ (^2) 



For the calculation of the spatial matrix elements we use the technique as described in 30 1 
to get 



{Fj'C'Mc' \'^l\FjCMc) 



2C 



KM" 



_d_ 

Ml 



j'CM" 



_d_ 



F. 



jCM" 



M",ti=± ) 



where drint = $,f$,2d^id$,2dt, t = ^i-^V The underlines mean that the space points take the 
value 

= 0, ^ly = 0, = 6x = ^2^1^, i2y = 0, 6. = (B4) 



the orbital angular momentum is given by I = x tti, and 
with the derivatives taken at the space point of (lB4p . 



(B5) 
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Appendix C: Calculational Details for V{NN A^A) Potential 



In calculating matrix elements of the transition potentials between antisymmetric ba- 
sis states one may omit the factor (1 - (12))/v/(2(l + B'^) from Eq. by usine; the 
substitutions 



1-(12] 



v/2(l + 



-.V 



NA,NN 



^3 - V2Vi2iNN ^ NA), 



Each basis state in fl23l) is the sum of several terms of the form 

J^^iTb'^) I ^^^■^ ® ^^'''^^ ^''''^ ^^^^^ ' ^^^''^''^ ^'^''^ ™^ 

and this is also true for the basis of the NNN part, but with B = B' = 1, (si, 52,53) 
= (1/2,1/2,1/2), and (1 - (12)) / [2(1 + i?'^)] 

elements of the operator V12 we need 



1. Therefore for the matrix 



1 

z 



A-ir''^''^'^[^llll]{^^^^^^^^ (CI) 

In the spatial matrix elements entering here we note that the NNN component and the 
NNA component of the wave function are given in terms of the Jacobi vectors of the same 
form (1221) but with different mass numbers. To perform the integration one needs to express 
one set of the Jacobi vectors in terms of the other via 



^1(^1,^2) 
^2(^1,^2) 



^1 



' A'2A[{Ai + A2 
{A[ + A',)A^A2 

iA[A2 - A'^A,) 



, iA[ + A',)A',A , 



(C2) 



^/A'{A[ + ^'2)^1^2(^1 + A2) ^A'{Ai + A2)As 
where A' = A'-^ + A2 + ^3, A = Ai + A2 + A3. The integration is done as 

{Fyc\\Vo{mB)\\F^c) = 87r\2C+l)-'/^ j rfr,„i 5^£>^(|i, |2)K,(mB)^kM(Ci, ^2)- (C3) 
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In addition we need ( using J\f^ — 4(7r/5)^/^l2M(n)) 



SS' J' 

x(((s;4)^'4)5'||E(^)||((siS2)5.3)5), 



{Fya\\V2{mB)M^'^\F,c) 



(C4) 



and 



M 



{{{s',s',)S's',)S'\\i:^'^\\{{s,s,)Ss,)S) 



-1) 



30/3S'SS'S 



S' S2 
S S' s' 



s'l Si 1 

s'2 S2 1 } (s;ikilki)(4lk2||s2)(C5) 

S' S 2 



Appendix D: Tj^ Multipoles of One-Body Currents Relating A 

For the magnetic multipoles one has 

i 

we have 



rpj, spin 
AAJrn 




(D3) 
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-t A A ^^yi- J 



1 ^l 1 



■AAjmV / r/iA V2 ■ 2 

The quantity d'^ is defined by the relationship —id'^ = p'^, and 



5e 



2,fi 



For the electric multipoles 



rpl \ ^ rry-ii, Spill / -N , ry-iL COnv / • \ n 



(D4) 



(D5) 



where / = j ± 1. One obtains 



-'-ANJm \^) 



mA 2 



2j + l 



(D6) 



ry-ijltl, spin 

^ AA,jm 



0-,; 



(D7) 



rpj±l,conv 



jm 



m 



For the combined r- and spin space operator we only need the reduced matrix element 
for the calculation of the response function, as seen in The reduced matrix element is 

given by 

{{Ffc ® {{s[s'^) S's'^)S')j,\\ {Oi ® Os)j \\{F,c ® {{siS2) Sss)S)j) 
C C I 
S' S s 



J' J J < 



} {F,,a\\Oimc){{{s[s',) S's',)S'\m{{s^S2) Sss)S), (D9) 



J' J J 



where 



{F^'C'\\Oi\\Fjc) = 87r^i2C' + l)-''/^ j dr.^tj^^- 

Mm 



and Oj is a function of relative coordinate and momentum 



1/2 



.2) 



M P3 - -jPcm 



A 



1/2 



_d_ 

W2 
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